The MiniBooNE anomaly, the decay Df — > n + v u and heavy sterile neutrino 
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It has been recently suggested that the anomalous excess of low-energy electron-like events 
observed by the MiniBooNE experiment, could be explained by the radiative decay of a heavy 
sterile neutrino Vh of the mass around 500 MeV with a muonic mixing strength in the range 
f/ M h| 2 ~ (1 — 4) x 10 -3 . If such vu exists its admixtures in the decay Dt — > p^v^ would re- 
sult in the decay Df — > ^ + Vh with the branching fraction Br(D^ — + ^ + vh) — (1.2 — 5.5) x 10 -4 , 
which is in the experimentally accessible range. Interestingly, the existence of the Df — > fi + Vh de- 
cay at this level may also explain why the currently measured decay rate of —* y^v^ is slightly 
higher than the predicted one. This enhances motivation for a sensitive search for this decay mode 
and makes it interesting and complementary to neutrino experiments probing sterile-active neutrino 
mixing. Considering, as an example the CLEO-c experiment, we suggest to perform a search for the 
decay Df — > fi + ish with the analysis of existing data. The discrepancy between the measurements 
and theoretical description of the decay — > t + v t is also discussed in brief. 

PACS numbers: 14.80.-j, 12.20.Fv, 13.20.Cz 



I. INTRODUCTION 

The MiniBooNE collaboration, which studies the in- 
teractions of neutrinos from the ir + decays in flight at 
FNAL, has observed an excess of low energy electron- 
like events in the energy distribution of charge-current 
quasi-elastic electron neutrino events [1]. This anomaly 
has been recently further confirmed with larger statistics 
[2]. As the collaboration has not yet clarified the origin 
of the excess, several models involving new physics were 
considered to explain the discrepancy, see e.g. [2] and 
references therein. 

It is well known, that the neutrino weak flavor eigen- 
states (24, u^, v r ,...) need not coincide with the mass 
eigenstates (y\, v%, V3, ^4...), but would, in general, be 
related through a unitary transformation. Such a gener- 
alized mixing: 

v l=/l Uu v i] I = e,M,T,..., % = 1,2,3,4, ... (1) 

i 

results in neutrino oscillations when the mass differences 
are small, and in decays of heavy neutrinos when the 
mass differences are large. In the recent work [3] it 
has been shown that the MiniBooNe excess could be ex- 
plained by the production of a sterile neutrinos, i^'s, of 
the mass around ~ 500 MeV, which, being created by 
mixing in neutral-current interactions, decay (dom- 
inantly) into photons and light neutrinos in the Mini- 
BooNE detector target. Such kind of Vh could arise in 
many interesting extensions of the Standard Model (SM) , 
such as GUTs, Superstring inspired models, Left-Right 
Symmetric models, and others. It can decay radiatively 
into z/7, if e.g. there is a non-zero transition magnetic 



* Sergei. Gninenko@ccrn.ch 
f Dmitry. Gorbunov@cern.ch 



10 





Excluded 




Allowed 


Excluded 








MiniBooNE 






favored 














i i i i i 



300 400 500 600 700 

m v . MeV 



FIG. 1: The shaded area is the experimentally allowed region 
of the mixing strength |c7 Mh | 2 calculated for jj, tr = 10 _9 /is 
in the model of Ref. [3]. The rectangular area represents 
the region in parameter space favorable for the suggested in 
Ref. [3] explanation of the MiniBooNE anomaly. 

moment (fJ-tr) between the Vh and active neutrino v [4]. 
The required mixing strength 

IL^I 2 ~ (1-4) x TO" 3 (2) 

was found to be consistent with existing experimental 
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data for \i tr ~ (1 — 6) x 10~ 9 /ib (here /j,b is the Bohr 
magneton) [3]. For illustration, experimentally allowed 
region of the mixing strength IL^I 2 in the vt mass range 
around 500 MeV is shown in Fig. 1 for [i tr = 10~ 9 yU-B to- 
gether with the parameter region favorable for the expla- 
nation of the MiniBooNE anomaly. It worth to mention 
that the model [3] is also consistent with the absence of a 
significant low-energy excess in MiniBooNc antincutrino 
data [5]. 

In this letter we put forward an idea that sterile-active 
neutrino mixing in the allowed range shown in Fig. 1 
could be tested by searching for the admixtures of 
in the decay Df — > /i + z^. In addition we point out 
that the present discrepancy of about 3cr between the 
measured and predicted decay rate of Df — > could 
be explained by the unrecognized contribution from the 
decay Df — > ^ + Vh- 

II. THE DECAY Df -> /i+i/ M AND HEAVY 
NEUTRINO 

If the Vh exists, it could be a component of v^, and as 
follows from Eq. (1), would be produced by any source of 
according to the mixing |C/ M h| 2 and kinematic con- 
straints. In particular, Vh could be produced in any 
lcptonic and scmilcptonic decays of sufficiently heavy 
mesons and baryons. For the interesting mass range 
m Vh ~ 400 — 600 MeV the most promising process is 
the leptonic decay Df — > ^ + v^. 

In the SM, D s meson decays leptonically via annihila- 
tion of the c and s quarks through a virtual W + . The 
decay rate of this process is given by 

nni - = - Jj^j \v ca \\ (3) 

where the M D + is the Df meson mass, mi is the mass 
of the charged lepton, f D + is the decay constant, Gf 
is the Fermi constant, and V cs is a Cabibbo-Kobayashi- 
Maskawa matrix element which value equals 0.97334 [6]. 
The decay rate (3) is suppressed by the lepton mass 
squared, since the very leptonic decay is due to chirality- 
flip. 

The mixing between the sterile neutrino and muon 
neutrino results in the decay Df — > fJ, + Vh, as illus- 
trated in Fig. 2. For the interesting mass interval m lJh ~ 
400 — 600 MeV the chirality-flip is mostly due to sterile 
neutrino mass which results in 

T(Df - »+v h ) « r(D+ - ^^)\U, h \ 2 (^) 2 . (4) 

Using Eq. (2) and taking into account the most pre- 
cise determinantion of the Df — > l^v^ branching ra- 
tio Br(L>+ -> ^+z/ M ) = (0.565 ±0.045 ±0.017)% [7], we 
find that the branching fraction of Df — > j-i + fh is in the 
experimentally accessible range: 

Br(^ ^ ^) « (1.2-5.5) x 10- (^^) 2 . (5) 




FIG. 2: Schematic illustration of the decay Df — > fi + Vh- 



III. DIRECT EXPERIMENTAL SEARCH FOR 
THE DECAY Df -> n+v h 

Consider now, as an example, the CLEO-c experiment, 
where the search for the decay Df — > i-i + v h could be 
performed. In this experiment several of the most pre- 
cise measurements of properties of Df mesons have been 
performed by using the CLEO-c detector at CESR [8]. 
Recently, the CLEO collaboration studying the process 
e + e~ — > DjD* + ,D*~Df has reported on measurements 
of the decay constant f D + of D s mesons to a precision of 
a few % [7, 9], see also [10]. 

The detector is well equipped to identify and measure 
the momenta/energy and directions of charged particles 
and photons. The experiment was performed at a centre- 
of-mass energy of 4170 MeV, where the cross section of a 
charmed meson pair production is relatively large. This 
allowed to fully reconstruct the D~ as a 'tag' and study 
the leptonic decay properties of the other through the 
decay chains 

e+e- -> D-D*+ - D~ 1 Df -> D~ sll x + v 
e+e" - D*-Df -> 7 D~Df -> -/>./'/' 

The decays Df — ► were identified by selecting the 

events with a single missing massless neutrino, for which 
the missing mass-squared, MM 2 , evaluated by taking 
into account the reconstructed n + , D~ and 7 should 
peak at zero. The MM 2 is calculated as 

MM 2 = {Ecu -E^-E 1 - E Ds ) 2 

- (PCM -Pn-Pl -PD S ) 2 

where Ecu and pcm are the centre-of-mass energy 
and 3- momentum, Ejj e and po s are the energy and 3- 
momentum of the fully reconstructed Dj tag, Ej and p 7 
are the energy and 3-momentum of the photon and E^ 
and Pn are the energy and 3-momentum of the muon [7] . 

Similarly to this approach, the basic idea of probing 
the model under discussion is to search for a peak corre- 
sponding to the value m 2 h in the MM 2 distribution. It 
should be calculated taking into account measured prop- 
erties of observed jj,, Df and photon from the decay 
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chains: 



these experiments one arrives at 



e + e 



e + e 



d;d;+ - Dj>yD+ - Djw+Vh - 

£>;-£)+ - 7 D a -£>+ - -tD-^+Vh - 

7^7 M + 7/^ 



(8) 



where 7^ denotes a photon from the dominant decay 
mode 1//, — > 7^ of sterile neutrino. In the largest part of 
the (|i/ M ft| 2 , ^i r ) parameter space favored by MiniBooNE, 
the Vh is expected to be a short-lived particle with the 
lifetime less than 10~ 9 s [3]. Then, its decay length is 
significantly less than the radius of the CLEO detector 
(95 cm), and most of the — ► 7^ decays would occur in- 
side the CLEO-c detector fiducial volume in the vicinity 
of the primary vertex. 

The experimental signature of the decay — > jis is a 
peak in the mass range 0.16 — 0.36 GeV 2 of the distri- 
bution of (7). Using Eqs.(2),(4) and the total number of 
(235.5 ± 13.8) M + t-V events observed by CLEO-c with 600 

pb _1 [7], (5 — 22) x ( 50 "^ cV ) 2 events are expected to be 
found at the peak. To obtain the correct MM 2 value, 
the photon from the decay — > 7^ should not be used 
in calculations of Eq. (7). For the energy greater than 
threshold, E lh > 300 MeV, the 7/, could be identified 
as an extra photon in the event candidate for the decay 
chains (6). 

Finally, note that a search for the decay mode Vh — ► /X7r 
is also of a special interest. Although this decay is sub- 
dominant, for the mixing as large as in Eq. (2) its branch- 
ing fraction could be of the order of few % [3] , and a few 
events could be observed in the CLEO-c experiment for 
the 600 pb _1 of data. The experimental signature of the 
event Vh — ► M 71 " would be two charged tracks originated 
from a common vertex displaced from the primary ver- 
tex. Since there is no neutrino in the final state, it is 
possible to reconstruct the invariant mass of the heavy 
sterile neutrino, that would manifest itself as a peak in 
the range 0.16-0.36 GeV 2 of the invariant mass squared. 
An observation of a few /X7r-events with the same invari- 
ant mass would provide an excellent cross-check of the 
model. 



IV. THE Df -> n + Vp,T + v T DECAYS PUZZLE 

Interestingly, the above discussions might be relevant 
to the discrepancy of about 3cr between the measured and 
predicted rates of Df — ► fj, + v^, see e.g. Refs. [7, 9] and 
discussion therein. Presently, in spite of the substantial 
theoretical and experimental efforts, the decay rates of 
Df — > p^v^ measured by the CLEO-c [7], Belle [11], and 
BABAR [12] experiments are found to be slightly higher 
than the predicted one, most accurately calculated in the 
framework of lattice QCD [13]. Taking into account the 
most precise results on the Df — ► decay rate from 



r e *P(£>+ -> /i+ty) - V th {Df -> 

T™P(Df -f ' (9) 

= 0.166 ±0.060 

where T ex P(Df -> fi+v^) and V th (Df -> n+v^) are 
the average measured and predicted values for the 
Df — > ^. + Vfji decay rate, respectively, with the statistical 
and systematical uncertainties combined in quadrature. 
Thus, the ratio of Eq. (9) differs from zero by about 2.8a 
standard deviations. 

Various models of new physics giving additional contri- 
bution to the rate of Df — > p^v^ have been investigated 
in order to resolve the discrepancy, see e.g. [14, 15]. We 
propose here that the reason of why the experimental 
rate of Df — > /i + iy is higher than the theoretical ex- 
pectations may be due the contribution from the decay 
Df -> n+v h . 

Consider again, as an example, search for Df — > fi^v^ 
events in the CLEO experiment [7]. If the 7^ from the 
decay chain (8) is used, the calculated MM 2 should peak 
at zero regardless of wether or not the 7^ is produced in 
the direct Df decay. This is valid under assumption 
that the v from the Vh — > 71/ decay is light. Thus, the 
events (8) may be accepted and contribute to the number 
of the Df — > n + v^ signal events. Using Eqs. (3,4,9) 
one finds that in order to explain the discrepancy the 
branching fraction of decay mode to sterile neutrinos of 
masses m Vh ~ 400 — 600 MeV should be within the range 



(9.85 ±4.16) x 10~ 4 



(10) 



where factor k = ^^i^ , Vh \ is the ratio of the overall 



P(D 



probabilities for the events (8) and (6) to pass selection 
criteria in the analysis of Df — > fJ^v^ events in [7]. In 
this search, for the selection of Df — > [i + v^ candidates it 
was required that there should be no additional photon, 
not associated with the tag, detected in the ECAL with 
energy greater than 300 MeV (photon veto) [7]. The 
fraction of 7^ from (8) that would pass this veto cut is 
roughly estimated to be k ~ 40%. This results in 



Br(£>+ -» ii+ vh) = (24.6 ± 10.4) 



x 10 



■ -4 



(11) 



The corresponding mixing strength is (c.f. (5) and (2)) 



\U^ h \ 2 = (25.2 ±10.7) x 10- 3 ^ 



1(K,M( ' V ) 2 . (12) 



The explanation of the MiniBooNE anomaly implies 
somewhat smaller contribution to the decay Df — ► 
However, the regions for the obtained mixing strength 
of Eq. (12) and branching fraction of Eq. (11) overlap 
at the level of less than 2 o with the ranges of Eqs. (2) 
and (5), respectively, required to explain the MiniBooNE 
anomaly, see also Fig. 1. This enhances motivation for a 
sensitive search of the decay Df — ► /i + z^- 



4 



In close analogy with the case of Df, the mixing (2) 
can also be probed with study of the lcptonic decay rates 
of the _D + -meson. In particular, our model implies the 
quite similar to (4) contribution to the total muonic decay 
rate of D + — > ^ + u. However, the achieved accuracy in 
measurement of this decay rate [6] is worse than that for 
the D s -meson, hence the suggested new contribution is 
unrecognizable yet. 

It should be mentioned that the best precision in mea- 
surements of the decay rate of Df — > t + v t has aslo been 
achieved at CLEO [9], and the result is consistent with 
theoretical predictions. However from other measure- 
ments at CLEO [7], based on study of another decay 
mode of outgoing r-lepton, the obtained decay rate of 
Df — > t + v t is significantly higher than the predicted 
one, so that combined branching ratio deviates from the 
theoretical prediction at the level comparable to Eq. (9). 
Thus, one may speculate that this discrepancy is due to 
existence of an additional sterile neutrino coupled to the 
tau neutrino. Similar to above consideration we found 
that the mixing strength required to explain Df —> t + v t 
discrepancy should be 

\U Th \ 2 = 0.16 ±0.09 (13) 

for the mass range much below 190 MeV, where phase 
space suppression is negligible. Note, that this result is 
consistent with direct experimental limits [16, 17]. 

It is worth noting that heavy sterile neutrino of mass 
about 500 MeV and mixing to muon neutrino as strong 
as (2), (12) can be searched in other decays of charmed 
hadrons, beauty hadrons and T-leptons. Pure leptonic 
modes are more promising, since heavy neutrino contri- 
bution is enhanced with respect to that of the SM by a 



squared mass ratio of sterile neutrinos and charged lep- 
ton. This is not the case for other decay modes which 
branching ratios are roughly the same as corresponding 
decay to active muon neutrino multiplied by \U^h\ 2 , and 
somewhat suppressed due to reduced phase space vol- 
ume. By making use of general formulae for meson and 
r-lepton decays to sterile neutrinos presented in Ref. [18] 
and similar formulae for baryons [19] one can obtain more 
accurate estimates. 



V. SUMMARY 

To summarize, we show that the recently suggested 
explanation [3] of the MiniBooNE anomaly [1,2] can be 
probed by studying leptonic decays of charmed mesons. 
This study can be undertaken with already collected 
data. We speculate that the present 2.8 a discrepancy 
between the measured and predicted rate of Df — > /U+z^ 
can be a hint at the presence of heavy sterile neutrino 
suggested in Ref. [3]. 
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